Probing quasi-particle states in strongly interacting atomic gases by 
momentum-resolved Raman photoemission spectroscopy 



o 
o 



C/2 



o 



00 

p 

in 
o 

O 



X 



Tung-Lam Dao,^ lacopo Carusotto,'^''^ and Antoine Georges^^** 

^Centre de Physique Theorique, Ecole Poly technique, CNRS, 91128 Palaiseau Cedex, France 
'"^CNR-INFM BEC Center and universita di Trento, 38050 Povo, Italy 
^Institute for Quantum Electronics, ETH Ziirich, 8093 Zurich, Switzerland 
^College de France, 11 place Marcelm Berthelot, 75005 Paris, France 
(Dated: May 6, 2009) 

We investigate a momentum-resolved Raman spectroscopy technique which is able to probe the 
one-body spectral function and the quasi-particle states of a gas of strongly interacting ultracold 
atoms. This technique is inspired by Angle-Resolved Photo-Emission Spectroscopy, a powerful ex- 
perimental probe of electronic states in solid-state systems. Quantitative examples of experimentally 
accessible spectra are given for the most significant regimes along the BEC-BCS crossover. When 
the theory is specialized to RF spectroscopy, agreement is found with recent experimental data. The 
main advantages of this Raman spectroscopy over existing techniques are pointed out. 

PACS numbers: 67.85.-d, 03.75.Ss, 42.65.Dr, 71.27.+a 



I. INTRODUCTION 

The recent advances in the cooling, trapping, and ma- 
nipulation of ultra-cold atomic gases using mostly op- 
tical beams have given birth to the new field of "con- 
densed matter physics with light and atoms" . Key issues 
of the physics of strongly correlated quantum systems 
can now be addressed in this new context from a com- 
pletely different perspective. Remarkable milestones in 
this respect have been the observation of the superfluid 
to Mott-insulator transition in a system of bosons in an 
optical lattice [l| , the direct imaging of the Fermi surface 
in a degenerate Fermi gas ^] , the demonstration of su- 
perfluidity in an interacting Fermi gas [1, 0, S H 0] ' ^^^"^ 
the recent observation [1, Q of an incompressible Mott 
insulating regime of fermionic atoms [l^, EH trapped in 
the periodic potential of optical lattices. 

Central objects in the theoretical and experimental 
study of quantum many-body systems are the low-energy 
excited states, and the possibihty of describing those 
states in terms of quasi-particle excitations having a long 
lifetime and a well-defined dispersion relation for the 
excitation energy as a function of momentum p^ . [isj . 
There is in fact abundant experimental and theoretical 
evidence that quasi-particle states can be highly uncon- 
ventional in strongly correlated systems, and significantly 
depart from the ones predicted by Landau Fermi-liquid 
theory [ij, 113 • An understanding of quasi-particle states 
is therefore an essential step in the direction of building 
a complete description of the peculiar electronic, mag- 
netic, and optical properties that have been observed in 
a variety of strongly correlated materials. 

A commonly used probe of the electronic states in solid 
materials is the so-called angle-resolved photoemission 
spectroscopy (ARPES) [l^. This technique consists in 
measuring the energy and momentum distribution of the 
electrons that are emitted from the solid when it is ex- 
posed to a beam of energetic photons. In the simplest 



approximation, the distribution of the emitted electrons 
is in fact proportional to the one-body spectral function 
of the electrons in the solid. 

Inspired by the success of ARPES in solid-state sys- 
tems, we recently proposed momentum-resolved stimu- 
lated Raman spectroscopy as a probe of quasi-particle 
states in strongly correlated atomic Fermi gases Ell ■ 
In a stimulated Raman process, atoms are transferred 
from the gas into a different internal state by absorbing 
a photon from a laser beam and immediately reemitting 
it into another beam of different frequency and wavevec- 
tor. For a given intensity and duration of the Raman 
pulses, the momentum distribution of extracted atoms 
is measured as a function of the wavevector and fre- 
quency of the Raman beams: analogously to electronic 
ARPES, this distribution is proportional to the one-body 
spectral function of the atoms in the gas. Alternative 
schemes that use Raman processes to probe the micro- 
scopic properties of atomic gases have been proposed 
in [ii, 113, iH, m m, m [ill. Experimental applica- 
tions of the related Bragg scattering technique have been 
reported in [H, 113, [H, . 

This momentum-resolved stimulated Raman technique 
is a powerful extension of the RF spectroscopy technique 
that has been recentl y u sed to study strongly interacting 
atomic Fermi gases jSOl . IsTI . [s^ . A pioneering experi- 
mental demonstration of the application of momentum- 
resolved RF spectroscopy to a strongly correlated atomic 
gas has recently appeared [11]: the energy dispersion 
of the occupied quasi-particle states has been measured 
at several points along the BEC-BCS crossover and has 
shown clear evidence of a pairing gap. Momentum- 
resolved RF spectroscopy in atomic Fermi gases has also 
been discussed recently in Ref. [s^ . 

In this paper we extend the discussion of Ref. 
and we give a comprehensive discussion of the promise 
of momentum-resolved stimulated Raman spectroscopy 
to investigate the microscopic properties of strongly in- 
teracting Fermi gases. On one hand, RF spectroscopy 



naturally arises as a special case of Raman spectroscopy 
with a vanishing transferred momentum: the theoreti- 
cally calculated spectra are in good agreement with the 
experimental observations of Ref.'ss'. On the other hand, 
Raman spectroscopy offers several significant advantages 
as compared to the RF one, e.g. spatial selectivity to 
eliminate inhomogeneous broadening due to the trap po- 
tential, tunability of the transferred momentum from be- 
low to well above the Fermi momentum, weaker sensitiv- 
ity to final-state interactions. 

The paper is organized as follows: in Sec|TTl we re- 
view the general theory of the Raman spectroscopy and 
we provide explicit formulas for the observable signal in 
terms of the one-body spectral function of the gas. In 
Sec lIIIi the advantage of the Raman spectroscopy over 
the RF one are illustrated on the simplest case of an 
ideal Fermi gas. In Sec lIVi we show how the general the- 
ory is able to reproduce the experimental observations 
of Ref 33 for a strongly interacting Fermi gas along the 
BEC-BCS crossover. The promise of an inverse Raman 
spectroscopy to investigate the dispersion of the empty 
quasi-particle states is discussed in SecEl Concluding 
remarks are given in Sec lVIl 

II. GENERAL THEORY OF THE RAMAN 
PHOTO-EMISSION TECHNIQUE 

A. The optical process 

In this section, we review the general theory of the 
stimulated Raman spectroscopy technique to probe the 
one-particle excitations of a many-atom gas that was 
originally proposed in Throughout the whole pa- 

per we will concentrate our attention on the case of a 
two-component mixture of fermionic atoms in two inter- 
nal states a and a', but the extension to a generic system 
of fermions or bosons is straightforward. Generally, the a 
and a' states are long-lived hyperfine components of the 
electronic ground state of the atom. Atoms are trans- 
ferred from the a state to a third hyperfine component 
/3 by a pair of laser beams in the Raman configuration 
schematically shown in Fig[TJ The laser fields are as- 
sumed to be far from resonance with the intermediate 
excited state 7 so that spontaneous emission events can 
be neglected. 

The atom-laser interaction can be described by the 
standard dipolar Hamiltonian 

Fdip = -d-[Ei(r,<)+E2(r,<)]. (I) 

where the laser fields are assumed to be classical and have 
the form 

Ei,2(r, t) = Ei,2(r) e*(ki.2r-..i.2t) _^ ^.c. (2) 

The envelopes Ei^2(r) are slowly varying on the scale of 
the carrier wavevectors ki.2 and account for the trans- 
verse profile of the laser beams. The corresponding 
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FIG. 1: Upper panel: geometrical configuration of the Ra- 
man laser beams. Lower panel: scheme of the atomic levels 
involved in the Raman process. 

Rabi frequencies are defined in terms of the electric 
dipole matrix elements of the optical transitions con- 
necting the a, (3 states to the common excited state 7 
as hni{r) = d-yQEi(r) and hn2{r) = d-y^E2(r). 

Once we perform the rotating wave approximation and 
we eliminate the intermediate excited state 7 under the 
condition fli,2, F-^ <C A (F^ is here the natural linewidth 
of the excited 7 state), the light-matter interaction can 
be reduced to the simple form: 

V{t) = Ve-''^' + h.c. (3) 

in which V — J dr fle{r,t) ipj^^-ipar e^^''^ ■ The basic pro- 
cess described by the Hamiltonian ^ consists of the 
transfer of atoms from the a to the 7 state by conser- 
vative Raman processes. ile{r,t) = r2i(r, i) 02(r, i)/A is 
the Rabi frequency of the effective Raman coupling, q = 
ki — k2 is the transferred momentum, and uj = lu i — UJ2 
is the transferred energy. 

Note that RF spectroscopy with a single electromag- 
netic field at frequency uj is described by a Hamilto- 
nian of exactly the same form Q with a Rabi frequency 
rie(r) = d^aE(r) and a vanishing transferred momentum 
q = 0. All the discussion that follows is then straighfor- 
wardly extended to the case of RF spectroscopy by simply 
setting q = 0. 

B. The Raman emission rate 

We first calculate the Raman transfer rate in the sim- 
plest case in which atoms in the final state (3 do not in- 
teract with the atoms left behind in the a, a' states and 
furthermore do not feel any trapping potential. Under 
these assumptions, their momentum is conserved while 
they freely propagate in space and one can envisage to 
perform a momentum-resolved measurement along the 
lines of the RF experiment of [3^ . 
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In terms of the many-body states of the gas, the Ra- 
man process then consists in the excitation of the ini- 
tial N body state lipf) to an excited — 1 body state 
and the simultaneous out-coupling of one atom 
into the /3 state with momentum k. The rate Rq(k,uj) 
for this process depends on the wavevector q = ki — k2 
and frequency w = 0i — 02 that are transferred by the 
Raman beams to the atoms. The total Raman trans- 
fer rate i?q (w) is then obtained by integrating over final 
momenta k as follows: 



(4) 



Direct application of the Fermi golden rule in the grand- 
canonical ensemble leads to 



5{Ef+ei,p-E,-huj). (5) 



The initial state is assumed to be at thermal equilib- 
rium at a temperature T with a chemical potential ^; 
Z — ^■exp{—Ei/kBT) is the corresponding Grand- 
Canonical partition function. The sums over i and / 
refers to all the many-body states of the system. The en- 
ergy Ei — Ei — fjiNi is rescaled by the number of particles. 
The matrix element has the form: 



rf^r(k/3|,/.j3jO)((/./|V'ar|(/'.)f^e(r) ^ 



(6) 



Using the definition of the real space spectral function at 
finite temperature 

A(v, r'; oj') = 

X {<f>^\^l,,\<t'f){^f\^o.r\<t>^)S{hl0' + Ef - E,) (7) 

the Raman transfer rate can be rewritten as 

i?q(k,cj) = ^ Jdh J^r' n,{r) r!:(r')e'(i-'^)('^-""') 
X nF(ek/3 - ftw - /i) A(r,r';ek/3 - ftw - m)- (8) 

Here, the Fermi factor has the usual definition npiE) = 
l/[l + exp{E/kBT)]. 

C. Spatially homogeneous geometry 

In the simplest case of a spatially homogeneous system 
with no trapping potential and spatially uniform Raman 



beams r2e(r) = ^e, the Raman rate per unit volume can 
be rewritten in a simplified form: 

dRqiKuj) 27r 2 / N 

— ^— = — \ile\ npiekfi -hhJ-^i) 

dV n 

X A(k-q,ek/3-ftw-M)- (9) 
that only involves the momentum space spectral function 



E 



^-Ei/ksT g-Ef/kBT 



X\{^f\lb^k\^,)\' Sihw + Ef ~E,). (10) 



In the simplest case of a non-interacting gas at zero tem- 
perature r = 0, the spectral function pH)) has the form 



v4(k, uj) ~ 6{fujj + fi — Ska) 



(11) 



The energy and momentum conservation condition 
that underlie ([9|) have a simple physical interpretation. 
Initially, the many-body system is in the state i of en- 
ergy Ej^ of momentum kf^; the laser beams consist of 
Ni^2 photons in the two beams of respectively frequen- 
cies 01^2 and momenta ki 2. At the end of the Raman 
process, the many-body system has an energy E^~^ and 
a momentum k^~^, while the laser beams contain re- 
spectively A^i — 1 and 7V2 + 1 photons. By energy and 
momentum conservation, one therefore has that 

-I- NihjjJi + N2huj2 = Ef-^ + {Ni - l)hui + 

+ (^2 + 1)^2 + ek/3 (12) 

kf + iViki -f 7V2k2 = kf -1 -f (TVi - l)ki + 

+ (iV2 + l)k2 + k. (13) 



which reduces to: 



N ,N-1 



= ek0- hu! 
= k~q 



(14) 



As expected, these conditions perfectly correspond to the 
frequency and momentum values at which the spectral 
function in © is evaluated. 



D. Local density approximation 

In view of concrete applications, it is useful to derive 
approximate formulas that can accurately describe the 
case of spatially selective Raman processes in trapped 
systems. Spatial selectivity is in fact a key advantage 
of Raman techniques over RF spectroscopy, as it allows 
to avoid inhomogeneous broadening effects by restrict- 
ing the Raman out-coupling process to a limited region 
of space where the system can be seen as almost uni- 
form. This possibility is even more important if the sys- 
tem presents several different phases with macroscopi- 
cally different properties. 
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The simplest way to include the trapping potential 
Vq,(R) acting on the a atoms and/or of the spatial pro- 
file of the Raman beams is to perform the so-called Lo- 
cal Density Approximation (LDA). This approximation 
is generally accurate as long as the characteristic length 
on which the properties of the system vary is much larger 
than all microscopic scales of the system. 

Under this approximation, the Raman rate can be 
written as an integral over different regions of space, 

X A(k - q, ek/3 - - Mr) (15) 

the contribution of each spatial region is approximated 
by the prediction ^ for a homogeneous system with an 
effective chemical potential /ir — /i— Va(R) that includes 
the effect of the trap. Here, A(k, hw; mr) is the spectral 
function (jTU]) with the local chemical potential /xr- 

By pushing further the LDA approximation, one can 
obtain a simple formula for the Raman rate also in the 
presence of an external potential Va(R) acting on the 
atoms in the f3 state, 

i?q(k,w) = ^ J d^Il\n,\^ npie^p - huj - fi) 

X ^(k- q,ek/3 - M;Mr)- (16) 



The only difference with (jisp is that the energy dispersion 
of P atoms now depends on position according to 



-R 



ek/3 + l^/3(R) 



(17) 



It is important to stress that the validity of (|16p requires 
more stringent conditions on the measurement process 
than a standard LDA. First, the out-coupled (3 atoms 
must not significantly move during the Raman process 
nor be accelerated. Second, the trap potential for f3 
atoms should be switched off as soon as possible after 
the Raman process in order to minimize the distortion 
of the k-space distribution of the out-coupled atoms in- 
duced by the evolution in the trapping potential. 

Equation (fTB|) is the key result of this section and will 
be the base of the calculations that we shall present in 
the following Sections. 



E. Localized probe 

For the sake of simplicity, we restrict our attention to 
the case of a pair of Raman beams at an angle with 
transverse Gaussian envelopes (FiglT]), 

Ei(r) - E°e-"'/2a'-^'/2a' (18) 
E2(r) = ^2 6^^^/^'^^"^'^'''"'^^^'^°''"*^^/^'^^. (19) 

The Raman coupling amplitude has therefore a Gaussian 
shape localized around r = 



with a peak amplitude 



and a Gaussian matrix 



A 



sin (f> sin (f> cos (f> 
sin (j) cos (j) 1 + cos^ 
2 



(21) 



(22) 



Inserting this formulas into the general expression for the 
Raman rate ([8]) and then changing to relative p = r — r' 
and center-of-mass variables R = (r -I- r')/2, one obtains 



X e'^"^ ^'"'npieufi - hui - p,) A(R, p; £k/3 - huj - p). 

(23) 

For the sake of completeness, it is interesting to briefly 
discuss the case of a spatially homogeneous system in 
the absence of any trapping potential that is probed by 
tightly focussed Raman beams. This formula can be 
further simplified by performing the Gaussian integrals, 
which leads to 



Rq(k,w) 



(20) 



fi,det[L] 

X nF(ek/3 — huj — p) A{k. — q', £k/3 — huj — p). (24) 

As compared to the result ([9]) for the spatially homoge- 
neous geometry, the momentum distribution of the Ra- 
man out-coupled atoms is now smeared out by the convo- 
lution with the Fourier transform of the Raman extrac- 
tion profile. As expected, this latter has a width of the 
order 1/cr. 

A faithful image of the spectral and momentum fea- 
tures of the system is then obtained simply by choos- 
ing a value of a much larger than all microscopic length 
scales of the system. This condition is compatible with 
the spatial selection of an almost homogeneous region as 
soon as the system is macroscopic, i.e. has a much larger 
size than all microscopic length scales. 

It is interesting to note that the Gaussian factor in 
p4)) tends to a delta- function in the cr — > oo limit, which 
perfectly recovers the results of Sec lII CI for the spatially 
homogeneous system. On the other hand, all information 
on the momentum distribution is lost in an extremely 
localized Raman probe of size a <^ kp^ which provides 
momentum- integrated images for a given w. 



III. NON-INTERACTING FERMIONS: RF VS. 
RAMAN 

Once we have established a general theory of Raman 
spectroscopy, it is important to validate its performances 
on some specific examples whose physics is well-known 
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and under control. In this Section, we shall start from 
the simplest case of a non-interacting degenerate Fermi 
gas. In particular, we shall discuss the effect of the trap- 
ping potential and we shall demonstrate the advantages 
of Raman spectroscopy over RF spectroscopic in order to 
obtain spatially resolved information on the microscopic 
properties of the system even in the presence of a trap. 

Throughout the whole section we shall use the pa- 
rameters of the Boulder experiment [s^: (i) a Fermi 
energy Ep — h x 9.4 kHz, (ii) a Fermi wave vector 
kp = 8.6 /xm~^, a total number of atoms = 3 x 10^, 
(iv) a total size of the cloud L = IdO/kp, and (v) a tem- 
perature T/Tp = 0.18. As the total size of the cloud 
L = 190/A:f is much larger than the only microscopic 
length scale kp^, the system can be safely considered in 
the macroscopic regime where the Local Density Approx- 
imation is accurate. 

In order the help the reader to extract the significant 
information on the dispersion of a-state quasi-particles, 
all plots from now on (except the left panel of Figd]) will 
be given as a function of = e^p — fku and k = k — q. 

A. RF spectroscopy 

Simulated plots for the RF signal are given in Fig[2] 
for different trapping configurations. As usual for RF 
spectroscopy, the RF field acts uniformly on the whole 
atomic cloud and the transferred momentum is q = 0, 
which implies k = k. 

(a) (b) (c) 




0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 
k/kp k/kp k/kp 

FIG. 2: Contour plot of the calculated RF intensity for a 
degenerate gas of non-interacting Fermions at T/Tf = 0.18. 
The color scale indicates the rate of transfer into a state of 
momentum k of the f3 state as a function of the (rescaled) 
RF frequency Es — ek/3 — i^. The left (a) panel refers to 
a spatially homogeneous system. The (b,c) panels have been 
calculated within the Local Density Approximation: The cen- 
tral (b) panel refers to the case where both a and /3 atoms 
experience the same harmonic trap potential. The right (c) 
panel refers to the case in which only the a-atoms are trapped. 
The RF field is spatially uniform over the whole system. 

The left panel Fig[2I^a) refers to the case of a spa- 
tially homogeneous system in the absence of any trap- 
ping potential. As the dispersion of the a and P states 
of non-interacting Fermions are exactly parallel, the 
momentum-resolved RF signal peaks at the same value 
tUres = Ek/J — Ska independently of k. Once we move to 



the (k, Es) variables considered in the plot, the peaks lies 
on the dispersion Ska of the a state. The peak intensity is 
independent of k and extends up to the edge of the Fermi 
sphere at kp. The abrupt edge at fc^^^ is here smoothened 
by the finite temperature value T/Tp — 0.18. 
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FIG. 3: Momentum-integrated RF signal for a trapped g 
in the central panel of Fig[2jb) , the trap potential is assumed 
to act in the same way on the two a, 13 atomic states. Differ- 
ently from (HI, the momentum integration is here performed 
along lines of fixed -Es. The calculated spectrum is in good 
agreement with the experimental observation by the Boulder 
group (Fig. 5b of Ref. 1331 ) 

The effect of trapping is included in the central and 
right panels, Figl2ljb,c). In the central panel, the a and 
(3 states are assumed to experience the same trapping po- 
tential as in the Boulder experiment. Under this condi- 
tion, the dispersions of the a and (3 states remain parallel 
at a fixed distance, so that it is still possible to extract 
the a state dispersion following the peak of the signal 
intensity in the (fc, Es) plane. In contrast to the previous 
case, the peak intensity is however strongly dependent on 
k: while the low-momentum states are filled at all posi- 
tions of the trap, only the center of the trap contributes 
in fact to the high- momentum states close to kp. This 
implies that the peak intensity strongly decreases with k 
and, in particular, it can be hardly visible in the region 
around kp. This can be a serious issue when one is to 
assess the effect of interactions, which is generally most 
pronounced in the vicinity of the Fermi edge. 

Note the qualitative agreement of the calculated sig- 
nal of Figl^Jb) with the experimental observations by the 
Boulder group shown in Fig. 3a of Ref. ^335 . The agree- 
ment is even more striking when we compare the theoret- 
ical plot for the k-integrated signal shown in Figl3] and 
the experimental result in Fig. 5b of Ref. [ssi- 

The effect of the trap is even more dramatic if we con- 
sider the case where only the a state feels the trap po- 
tential while the (3 atoms are untrapped. The distance 
between the and dispersions is now strongly de- 
pendent on the position R, which produces the sizable in- 
homogeneous broadening that is visible in Fig[2l^c): the 
lower boundary corresponds to the contribution of the 
central region of the trap, while the upper boundary cor- 
responds to the edges of the cloud. In this case, it ap- 
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pears to be difficult to extract useful information on the 
£ka dispersion from such a broadened RF signal. 



B. Raman spectroscopy 

The advantage of the Raman spectroscopy technique 
over the RF one is clearly visible in the simulated Ra- 
man signal that is plotted in Fig HI F or the parameters 
of the Boulder experiment of Ref. l33l a transverse size 
of the Raman beams a = \b/kp ~ 1.74 /im well satis- 
fies the inequalities L 3> cr ^ . In particular, we 
consider the Raman beams to be focussed onto the cen- 
tral region of the trap and to have an angle </> = 7r/2 
between them. Their wavelength is assumed to be in 
the optical range, e.g. A = 0.7 /im. The momentum 
transferred to the atoms during the Raman process is 
then |q| = 2 sin((/)/2)/A « 2/im"^, much smaller than 
the Fermi momentum q/kp « 0.2. 




0.5 1 1.5 

(k-q)/kp 



FIG. 4: Contour plot of the Raman rate for a trapped gas 
with the same parameters as in the Boulder experiment and 
a beam transverse size a — 15/ kp- Left panel; raw data in the 
{k,uj) plane. Right panel: Translated data in the (|k — q|,_Es) 
plane. 

Raw data for the momentum-resolved Raman rate as a 
function of w and k are plotted in the left panel Fig.|4fa). 
In order to facilitate physical understanding of the quasi- 
particle dispersion, the same data are plotted in the right 
panel as a function of the rescaled variables Es = ek/3 — ^ 
and k = k — q. Even though we are dealing with a 
trapped gas, the observed signal is now very similar to 
the one of the homogeneous system shown in Figl^fa): 
this proves the utility of the spatial selectivity of Raman 
spectroscopy in view of extracting information on the lo- 
cal microscopic properties of a trapped gas. The broad- 
ening due to the finite beam size a = Ib/kp is almost 
irrelevant on the scale of the graphs. 



IV. STRONGLY CORRELATED FERMIONS 

To better understand the advantage of Raman spec- 
troscopy over the RF technique, we now extend our 
analysis to a remarkable example of strongly interacting 
fermionic system. Inspired by the Boulder experiment of 
Reflssl. we consider (i) a Fermi gas in the unitary limit 



{kpas)~^ = at a temperature close to the superfluid 
critical temperature where the physics is dominated by a 
pseudogap, and (ii) a molecular gas in the BEC regime 
{kpas)~^ > 1. In both cases, we will show how the Ra- 
man spectroscopy is able to provide useful information 
on the quasi-particle spectrum of the system. 



A. Preformed pairs at unitary limit 

At a temperature just below the superfluid critical 
temperature T/Tc = 0.9, the superfluid fraction and the 
superfluid order parameter are still very small. Yet, a 
sizable pairing gap appears in the spectral function due 
to the presence of preformed pairs. Remarkably, such a 
pseudo-gap has a width of the order of the Fermi energy 
in the critical region and persists even at temperatures 
well above the critical point. 




FIG. 5: Left (a) panel: Pseudogap model for a Fermi gas in 
the unitary limit. The black line is the free particle disper- 
sion while the blue and red curves are the BCS quasi-particle 
branches. Central (b) and right (c) panels: contour plots 
of the RF rate for a spatially homogeneous system at low 
T/Tp = 0.05 and intermediate T/Tp = 0.4 temperatures. 

A reasonable picture of the spectral function in the 
presence of a pseudogap is obtained using the standard 
BCS theory with the pseudogap Aps replacing the stan- 
dard superfluid gap Asp, 

^(k, Lo) = ui 6{huj + ti-E+)+vl 5{huj + ^l-E-) (25) 

The quasiparticle energy dispersion E^ has the usual 
BCS form 



A2 



as well as the Bogoliubov coefficients and v^, 

£ka — M 



1 ± 



(26) 



(27) 
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The RF and Raman signals for a spatially homogeneous 
Fermi gas at unitarity are readily obtained by inserting 
the spectral function ([25|l into the general expression ([9]) . 

At very low temperatures ksT/ Aps ^ 1, only the sec- 
ond term contributes to (j25p : the density of unpaired 
atoms is in fact exponentially suppressed by the factor 
exp(-Aps/fcsr) < 1. As shown in FigEl^b), the RF rate 
then peaks on a single branch in the (fc, Es) plane corre- 
sponding to the lower BCS branch at E^. This branch is 
strong up to the edge of the Fermi sphere. The long tail 
at high momenta is due to the particle-hole mixing char- 
acteristic of BCS theory that is visible in the expression 
([27)1 of the Bogoliubov coefficient. 

When the temperature gets higher, a second branch 
appears that corresponds to the upper branch of the BCS 
dispersion at E^. The weight of this branch rapidly 
grows with temperature as exp^—Aps/ksT) <C 1. As 
one can see in Figl5l[c), the combined effect of the 
coefficient and of the Fermi factor makes the intensity 
of this branch to be concentrated in the region close to 
the pseudo-gap where most of the unpaired particles are 
found. 




0.5 1 1.5 
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FIG. 6: Photo-emission signal for a strongly interacting, 
trapped Fermi gase at unitarity {l/kpcis = 0) at a tempera- 
ture T = O.ISTf. The left panel shows the RF signal. The 
right panel shows the Raman signal for a spatially selective 
process with a = 15/kF- For a comparison of the left panel 
to experiment, see Fig. 3b of Refls^. 

A quantitative comparison to the experimental data 
shown in Fig. 3b of Ref. [s^ requires that we include in 
the theoretical model the effect of trapping. This is done 
in Figinija): both atomic states a and f3 are assumed 
to feel the same trap potential. As we are working at 
the unitary limit at temperatures much lower than the 
pseudo-gap energy, this latter can be approximated to be 
proportional to the local Fermi energy which gives: 



Ap,(r) = Ap,(0) 1- 



(28) 



The broadening of the line that is apparent at fc = 
results from the spatial variation of the minimum 
of the quasi-particle dispersion, Emin = ^J■ — Vair) — 
~ V^(r)]^ + A^. The second branch that around 
k/kp ~ 0.8 emerges from the Bogoliubov dispersion in 
the upwards direction is due to the reduced gap ampli- 




FIG. 7: Momentum-integrated RF signal for a trapped, 
strongly interacting Fermi gas at unitarity. Same system as 
in the left panel of Fig|6] As in Fig|3j momentum integra- 
tion has been performed along lines of constant Es- For a 
comparison to experiment, see Fig.Sc of Reflssl. 



tude in the outer part of the cloud and to the corre- 
sponding unpaired atoms. The qualitative agreement of 
theoretical data with the experiment is satisfactory. A 
more quantitative comparison between theory and ex- 
periment (Fig. 5c of Ref. [H) is successfully made on the 
momentum-integrated spectral density shown in Fig[71 

The advantage of the Raman spectroscopy technique 
over the RF one is apparent in Fig|S]^b), where the sim- 
ulated signal for a Raman probe localized in the cen- 
tral part of the trap is plotted. As expected, the in- 
homogeneous broadening effects disappear and the sig- 
nal closely follows the Bogoliubov branch. Furthermore, 
the intensity of the branch is everywhere determined 
by the Bogoliubov Vk coefficient and does not involve 
any spatial average. For the parameters considered here 
{T/Tp = 0.18, A ~ Ep), the intensity of the upper Bo- 
goliubov branch of unpaired atoms due to the thermal 
excitations is relatively small. 



B. Tightly bound molecules in BEC limit 



The same BCS theory that we have used in the previ- 
ous section to describe the unitary gas can be extended 
to the BEC side where the gas is constituted by tightly 
bound molecules. The spectral function has the form 
but the gap A and the chemical potential ^ ap- 
pearing in the quasi-particle branches ((26)l have to be 
calculated by means of the self-consistency relations of 
BCS theory 

In the simplest case far in the BEC limit kpCLs ^ 1 at 
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T ~ 0, one has: 



INVERSE RAMAN SPECTROSCOPY 



A 




(29) 
(30) 
(31) 



which corresponds to a pair of quasi-particle branches 



— Eh — £ka 



TO 



(32) 
(33) 



The upper branch corresponds to thermaUy excited un- 
paired atoms. As a consequence of the (repulsive) poten- 
tial of molecules, the branch is slightly blue-shifted with 
respect to the free atom dispersion. The lower branch 
corresponding to molecules is separated by the upper one 
by the binding energy Ei, of a molecule and has an in- 
verted dispersion with negative effective mass. As sug- 
gested in [33| , this peculiar fact has a transparent physi- 
cal interpretation: eliminating a a atom of wavevector k 
means that one has to break a molecule and eventually 
leaves the system with an unpaired a' atom at a momen- 
tum — k. Not only this requires the binding energy E(, 
of the molecule, but also the kinetic energy e_ka' of the 
unpaired atom. 




0.5 1 1.5 2 
(k-q)/kp 



FIG. 8: Photo-emission spectroscopy signal for (kpas)^ — 1 
on the BEG side at a temperature T/Tp = 0.5. Left (a) panel: 
RF signal for a homogeneous system. Gentral (b) panel: RF 
signal for a trapped system with the parameters of the ex- 
periment of Ref 33. Right (c) panel: Raman signal with a 
spatially selective probe a — Ib/kp focussed onto the cloud 
center. 



These two branches are visible in the plots of the RF 
and Raman signals for a molecular condensate at finite 
temperature that are shown in FiglHl As usual, the Ra- 
man signal with a spatially selective probe [panel (c)] is 
able to overcome the inhomogeneous broadening effects 
due to trapping that would otherwise disturb the RF 
signal [panel (b)] and almost recovers the RF signal of a 
homogeneous system [panel (a)]. 



All the discussion so far has been focussed on the 
case where the state P is initially empty: in this case, 
momentum-resolved photo-emission spectroscopy is able 
to provide detailed information only on those bands 
that are initially occupied. This fact is apparent in the 
npi^up — hu! — ii) factor multiplying the spectral function 
A in the equation for the Raman rate ©. 

In many cases of current experimental interest, it is 
however desirable to have experimental access to the dis- 
persion of empty bands as well. This is particularly inter- 
esting at very low temperatures, where only the "negative 
energy" states below the chemical potential are occupied, 
but much interesting physics is contained in the "positive 
energy" ones above the chemical potential. For instance, 
in the BCS-like case discussed in Sec lIVi observation of 
both E^ branches would provide unambiguous informa- 
tion on the amplitude of the superconducting gap. 

Information on the empty bands can be obtained if 
both direct a — > /3 and reverse /? — > a Raman processes 
can be induced by the same pair of laser beams. This 
is the case if some incoherent population of atoms is al- 
ready present in the (3 state at the beginning of the Ra- 
man experiment: the relative occupation of the initial 
and final states determines whether the direct or the re- 
verse Raman process will dominate. The resulting signal 
results from the difference of the two photo-emission and 
photo-absorption processes and is quantified by the rate 
of increase/decrease of the population of the k momen- 
tum state of the /3 level. Within the same approximation 
performed in the previous sections of the paper, the rate 
([9]) is easily generalized to include also reverse Raman 
processes: 



dV ~ 



fku 



- A{k - q, Si^f3 



-m) 

M/3)] 



(34) 



h 

The ideal Fermi gas in the f3 state is here assumed to 
have a thermal distribution at same temperature T as 
the gas in the a state, and a chemical potential Gen- 
eralization of to the case of an arbitrary occupation 
law is done by simply replacing the occupation factor 
npi^kfs — M/s)- In np{e-ki3 — ^p) = limit, equation 
(jM)) reduces to ®. 

In the q = case of RF spectroscopy, probing the 
upper BCS band requires that a sizable population is 
present in the (3 state up to high momenta. At low tem- 
peratures, this requires that the Fermi momentum kpp 
of the (3 state is larger than the one kp of the a state, i.e. 
that there are initially more atoms in the (3 state than 
in the a one. Examples of combined photo-emission and 
photo-absorption spectra of the q = case are shown 
in Figl5]for the BCS-pseudogap model of strongly corre- 
lated Fermions discussed in Sec lIVI and described by the 
spectral function (P5|) . Different panels refer to differ- 
ent values of the initial population of the (3 state (left to 
right) and to different temperatures (up to down). 
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0.5 1 1.5 
(k-q)/kp 



2 0.5 1 1.5 2 
(k-q)/kp 



FIG. 9: Photo-emission spectroscopy with q = for differ- 
ent temperatures T = O.OSTf (upper panels), T — OATf 
(lower panels) and for different initial populations of /3-atoms 
kp/s = (left panels), kp/s = 1.73fcF (right panels). In all pan- 
els, the purple background color corresponds to a vanishing 
photo-emission and photo-absorption intensity. The yellow 
color indicates the occurrence of direct, photo-emission pro- 
cesses from a into /3; the blue color indicates the occurrence 
of reverse, photo-absorption processes from /3 into a. 



The left panels correspond to an initially empty f3 
state: at T/Tp = 0.05 (upper panel), only the lower BCS 
branch at i?^ is visible as a positive, photo-emission sig- 
nal (yellow). At higher T/Tp = 0.4 (lower panel), also 
the upper BCS branch at appears in the spectrum 
again as a positive signal. The right panels correspond 
to a highly degenerate (3 state where the lowest (3 states 
have an almost unity filling kpp = 1.73kp. At low tem- 
peratures {T/Tp = 0.05), photo-emission from the E^^ 
is strongly suppressed by Pauli blocking, while the up- 
per band at E^ clearly appears as a negative, photo- 
absorption signal (blue). At higher temperatures, both 
bands are visible due to thermal broadening. 

The difficulty of having a high initial density of atoms 
in f3 state can be overcome by adopting a Raman scheme 
with a transferred wave vector q comparable to kp (see 
FigfTU)). In this case, the population in the almost full 
lowest (3 states can be effectively transferred into the 
empty a states in the most interesting region k hp 
around the superconducting gap. 

Simulated spectra for higher values of the transferred 
momentum q = O.bkp, kp (left to right) are shown in 
Fig-EHl At low temperature T/Tp = 0.05 (upper panel), 
the photo-emission processes from the lower E^ branch 
are suppressed by Pauli blocking in the momentum range 
[—kpf3 — q, kpp — q]. By contrast, the photo-absorption 
processes into the upper E^ branch clearly appear in 





-2-10 1 
(k-q)/kp 



2 -2 



-1 1 
(k-q)/kp 



FIG. 10: Photo-emission spectroscopy for different tempera- 
tures T = O.OSTf (upper panels), T = OATp (lower panels) 
and for initial population of /3-atoms kpp — 0.5kp. The trans- 
fer momentum are: q — 0.5kp (left panels) and q = kp (right 
panels) . Same color code as in Fig|9l 



this momentum range as a negative signal. At a higher 
temperature T/Tp = 0.4 (lower panel), both processes 
are broadened over the whole range of momentum by 
thermal effects. 

The comparison between q = 0.5kp (left panels) and 
q — kp (right panels) shows that the double branch struc- 
ture is most visible when the transferred momentum q is 
closest to the Fermi momentum kp: in this case, the 
value of the pseudo-gap amplitude of the BCS model can 
be easily extracted from the distance in frequency be- 
tween the two features. 



VI. CONCLUSIONS AND PERSPECTIVES 

In conclusion, we have proposed a momentum-resolved 
Raman spectroscopy technique that, in analogy to the 
Angle-Resolved Photo-Emission Spectroscopy of solid 
state physics is able to probe the one-body properties 
of an atomic gas. The power of this technique to mea- 
sure the dispersion of the filled quasi-particle states has 
been illustrated on a number of simple systems ranging 
from an ideal Fermi gas to a strongly correlated one and 
has been validated on existing experimental data. 

Several advantages over previous techniques (in par- 
ticular RF spectroscopy) have been pointed out: the use 
of focussed Raman beams enables one to eliminate the 
inhomogeneous broadening due to the trap potential by 
restricting the optical process to a limited spatial region. 
A large value of the transferred momentum can also help 
to purify the measured spectra by suppressing the effect 
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of interactions between the photo-emitted atoms and the 
rest of the cloud. 

A direct extension of the technique to the case where 
an incoherent population is already present in the final 
state of the Raman process is finally proposed: Raman 
processes take place in both directions and reverse Ra- 
man spectroscopy can be used to obtain information on 
the empty quasi-particle states lying above the chemical 
potential. 

Thanks to these remarkable features, we expect that 
momentum-resolved Raman spectroscopy will play an 
important role in the experimental characterization of 
the variety of novel quantum states of matter that can 
be obtained in systems of strongly interacting ultracold 
atoms, e.g. fermionic Mott insulator states, d-wave su- 
perconductors, anti-ferromagnetic states. 
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